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Anastasios Avgoustidis, 1 ']^] Edmund J. Copeland, 1 ']^] Adam Moss, and Dimitri Skliros 1 ^ 

1 School of Physics & Astronomy 
University of Nottingham, Nottingham, NG7 2RD, England 
(Dated: September 13, 2012) 

We present analytic expressions for the cosmic string unequal time correlator (UETC) in the 
context of the Unconnected Segment Model (USM). This eliminates the need to simulate the many 
thousands of network realisations needed to estimate the UETC numerically. With our approach 
we can compute the UETC very accurately, over all scales of interest, in the order of ~ 20 — 30 
seconds on a single CPU. Our formalism facilitates an efficient approach to computing Cosmic Mi- 
crowave Background (CMB) anisotropies for cosmic strings. Discretising the UETC and performing 
an eigen-decomposition to act as sources in the CAMB CMB code, the power spectrum can be calcu- 
lated by summing over a finite number of eigenmodes. A much smaller number of eigenmodes are 
required compared to the conventional approach of averaging power spectra over a finite number 
of realisations of the string network. With the additional efficiency and performance improvements 
offered by the OpenMP CAMB code, the time required to compute string power spectra is significantly 
reduced compared to the standard serial CMBACT code. The latter takes ~ 30 hours on a modern 
single threaded CPU for 2000 network realisations. Similar percent level accuracy can be achieved 
with our approach on a moderately threaded CPU (8 threads) in only ~ 15 minutes. If accuracy 
is only required at the 10 percent level and the CPU is more highly threaded, cosmic string power 
spectra are now possible in ~ 2 — 3 minutes. This makes exploration of the string parameter space 
now possible for Markov-Chain-Monte-Carlo (MCMC) analysis. 

PACS numbers: 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 

Results from the COsmic Background Explorer 
( COBE) [T] and Wilkinson Microwave Anisotropy Probe 
( WMAP) [2] satellites have provided strong support 
that close to scale invariant fluctuations, with an am- 
plitude of ~ 10 -5 , were generated by a period of accel- 
erated expansion (inflation) in the early Universe. This 
paradigm provides an extremely good fit to other cosmo- 
logical data, so by the principle of Occam's razor, having 
any additional source of fluctuations, such as topological 
defects, is deemed unnecessary. 

This line of reasoning, however, is only valid when 
inflation and defects are independent of each other. In 
many theories of the early Universe the two are inti- 
mately related. Cosmic strings, for example, are ubiqui- 
tous in models of supersymmetric [2H3 and brane infla- 
tion [BHS] ■ Inflation ends via a phase transition, and it is 
this symmetry breaking which results in strings HOP 
[12] . The amplitude of fluctuations created by strings 
can be comparable to the inflationary component, so it 
is crucial to account for them for a self-consistent theory. 

However, strings should not be viewed as "necessary 
evil" imposed on us by self-consistency. Instead, they 
offer an exciting window of opportunity for exploring 
High-Energy Physics, through their potential observa- 
tional effects. In particular, cosmic superstrings [T31H4] 
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- generically produced in string theoretic brane inflation 
models [7] [TTJ [HI [TS] - form complicated multi-tension 
networks with segments joining together in Y-shaped 
junctions, and theses characteristics depend quantita- 
tively on fundamental string theory parameters and 
compactification data. Thus, they offer one of the very 
few observational windows into string theory, and even 
the non-detection of strings can constrain fundamen- 
tal string theory parameters |16j . For the simplest net- 
works of "ordinary" field theory strings, non-detection 
results in a constraint on the energy scale of the string- 
producing phase transition, which marks the end of the 
inflationary phase. For recent reviews of the opportu- 
nities offered by cosmic strings for probing inflationary 
models and string theory, see [T71 [18] . 

In this paper we consider the observational signature 
of cosmic strings on the Cosmic Microwave Background 
(CMB). Strings create anisotropies by their motion 
along the line-of-sight (the Kaiser-Stebbins effect [T5] ) 
between ourselves and the last-scattering surface, as well 
as other related physics at the epoch of recombination. 
Since the formation of strings is a non-Gaussian pro- 
cess, there are associated non-Gaussian signatures in 
the CMB (see Ref. [20] for a recent review). In this 
work, however, we only consider the Gaussian compo- 
nent, that is the CMB power spectrum. 

The amplitude of string fluctuations are typically 
characterised by their dimensionless mass per unit 
length, G^i/c 2 . Current CMB power spectra data al- 
lows for G/i/c 2 <2-6x 10~ 7 at 2a confidence [HH23]. 
The range in this limit arises primarily due to two rea- 
sons: (1) the precise details of the evolution of the string 
network, i.e. its scaling properties; (2) modelling uncer- 
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tainties due to the dynamical range of simulations. A 
recent discussion of the various assumptions and prob- 
lems with each method can be found in Ref. [2"T] , 

No matter which method is used, the key intermedi- 
ate quantity of interest is the Unequal Time Correlator 
(UETC) of the string network [24]. There exists a sep- 
arate UETC for each of the scalar, vector and tensor 
contributions to the energy-momentum tensor. From 
this, one can then calculate CMB anisotropies using 
an Einstcin-Boltzmann code. The UETC can be esti- 
mated directly from simulations. An alternative method 
is to model the strings as an ensemble of uncorrelated 
straight segments, with a root-mean-square (RMS) ve- 
locity and correlation length chosen to match that of 
simulations. The latter approach is known as the Un- 
connected Segment Model (USM) [2SJ[2B], and by mod- 
elling strings as ID objects it avoids many of the issues 
associated with dynamical range. It is also easier to 
vary the network parameters within the context of the 
model and check their impact on observable quantities. 
The USM has disadvantages of its own however, which 
are discussed in Ref. [21]. 

Up to now the UETC of the USM has only been (in- 
directly) calculated numerically, by simulating an en- 
semble of source histories. The physics of the USM is 
relatively simple, however, and in this paper we show the 
UETC can be derived analytically. This is an important 
result for two reasons. The first is that, since the the 
UETC encodes all information of the string model, it is 
extremely useful to be able to evaluate this quickly and 
assess the impact of changing parameters, or to compare 
it with full numerical simulations of the string network. 
The second is that knowledge of the UETC allows the 
CMB power spectrum to be computed in a more effi- 
cient way, rather than averaging over source histories of 
the USM. 

Our formalism has particular relevance to the next 
generation of experiments which will measure the polar- 
isation signal of the CMB. Since strings generate vector 
and tensor fluctuations, they also create a B-mode po- 
larisation signal. In many theories of the early Universe 
which give rise to inflation and strings, the B-mode sig- 
nal from inflation can be much smaller than that from 
strings. This means the string B-mode signal could pro- 
vide a unique window to testing high energy physics at 
the ~ 10 16 GeV scale. In a future publication we will 
use the analytic form of the UETC to approximate, an- 
alytically, the string B-mode spectrum in terms of the 
underlying parameters. 

The outline of the paper is as follows. In section [H] 
we derive analytically the cosmic string UETC's. In sec- 



tion III we use these as inputs to an Einstcin-Boltzmann 
code to compute CMB temperature and polarisation 
power spectra. Our method dramatically improves com- 
putational time, making exploration of the string pa- 
rameter space now possible for Markov-Chain-Monte- 
Carlo (MCMC) analysis. Finally we conclude and out- 
line future directions of work. Throughout this paper 
we will work in relativistic units with c = 1. 



II. STRING UNEQUAL TIME CORRELATOR 

Perturbations generated by cosmic strings are qual- 
itatively different from those arising from inflationary 
vacuum fluctuations. While inflationary perturbations 
are passive, in that they are set as initial conditions 
(on super-horizon scales) and are then evolved through 
the "source free" linearised Einstcin-Boltzmann equa- 
tions, cosmic strings generate fluctuations which add a 
forcing term to these equations. Thus, strings actively 
source metric perturbations (on sub-horizon scales) at 
all times, and so their cosmological evolution is of key 
importance for predicting CMB anisotropies. Indeed, 
CMB spectra can be computed from the unequal time 
correlators (UETC's) of the network energy-momentum 
tensor components (in Fourier space), which are mani- 
festly functions of the time-dependent network parame- 
ters. 

In this section, we describe our analytic approach 
for computing the relevant UETC's for cosmic string 
networks. Starting from the (Fourier transform of 
the) energy-momentum tensor of a single string seg- 
ment, we review and update the Unconnected Segment 
Model (USM) ;25,:26] procedure for building the energy- 
momentum tensor components for a network of strings. 
Given the UETC's, these are then used as inputs to an 
Einstein-Boltzmann code to calculate CMB power spec- 
tra. This is a key step in improving the speed with which 
we can calculate spectra over existing string codes. 

The fact that strings are active sources allows them 
to produce significant vector perturbations, which nor- 
mally decay in passive models where they are not contin- 
uously sourced. In fact, it is the vector metric perturba- 
tion that provides the dominant contribution to the B- 
mode spectrum, which is of particular interest for future 
CMB polarisation experiments. Unlike inflation, where 
tensor modes are suppressed by a factor proportional 
to the slow-roll parameter, strings generally also pro- 
duce a non- negligible amount of tensor perturbations. 
We therefore need to consider all three contributions, 
scalar (S), vector (V) and tensor (T), and compute all 
the relevant UETC's. 



A. Evolving string networks 

Since cosmic strings are active seeds of perturbations, 
their cosmological evolution affects the CMB. In order 
to confront cosmic strings with CMB data, it is therefore 
of crucial importance to be able to accurately model the 
behaviour of string networks throughout cosmic history. 

Over the last couple of decades, significant progress 
has been made in this direction and a number of com- 
plementary techniques have been developed to describe 
the evolution of a wide range of string models (For re- 
views see Refs. QH ET] [28]). The dynamics of "or- 
dinary" string defects is fully described by non-linear 
field theories in cosmological backgrounds, which can be 
evolved numerically in lattice simulations 29 31 . For 
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cosmological applications, short scale effects are gener- 
ally neglected (even though there remains a question 
as to the validity of this, see for example Ref. [32]) 
and the networks can be approximated by a zero-width, 
Nambu-Goto string limit. With some additional input 
from field theory to describe the result of string inter- 
commutations [SHI EH] , this has allowed high-resolution 
simulations of Nambu-Goto strings in realistic cosmolo- 
gies |35H39| . Comparing numerical results with analytic 
approaches, it has also been realised that key features 
of the complicated network dynamics can be captured 
by surprisingly simple analytic models, describing the 
string network by only a handful of macroscopic vari- 
ables. 

In this approach, one starts with the assumption (jus- 
tified by simulations) that the string network has a 
random walk structure, so it can be characterised - 
on large scales - by a correlation length, L, quantify- 
ing the average length of string segments and the aver- 
age separation between theirj^] Similarly, string motion 
can be described by a root-mean-square (RMS) veloc- 
ity, v, of string segments. These two variables satisfy 
macroscopic equations of motion, derived directly from 
the Nambu-Goto action, after averaging over the world- 
sheet and introducing a phenomenological interaction 
term 



1 dL , 9n cv 
L^ = (1 + V) + 2L' 



dv 
Itt 



L 



- 2Hv 



(1) 
(2) 



where H is the Hubble function, H = a/a, a(t) is the 
scale factor and an overdot represents the derivative 
with respect to physical time t. The parameter c is 
a constant quantifying energy loss due to loop produc- 
tion and can be calibrated by comparison to simulations. 
Finally, the curvature parameter k indirectly encodes in- 
formation about the small-scale structure on strings and 
can be expressed in terms of v. For relativistic strings 
we have [43] 



k = 



2^2 fl-Sv 6 
l + 8v 6 



(3) 



Equations (T|3) form the Velocity-dependent One-Scale 
(VOS) model, and have been shown to be in remarkable 
agreement with Nambu-Goto and field theory simula- 
tions for the simplest of cosmic string models. Writing 
the correlation length in terms of a new variable e, 



and expressing the Hubble function in terms of the ex- 
pansion parameter (3, such that a(t) oc t&, one finds that 
this system has the scaling solution 



k(k 



4)9(1 - P) 
P{k + c) ' 



(5) 
(6) 



which is the generic attractor. In comoving units, the 
correlation length is then 



a 



(7) 



where r is the comoving horizon (or alternatively con- 
formal time). Since, in a given era (f3 = const, equal 
to 1/2 and 2/3 in the radiation and matter eras respec- 
tively), the physical horizon is dn = at = i/(l — jS), the 
relation between e and £ is 



£ = £(1-/3) 



(8) 



L = et. 



(4) 



For string networks of higher complication, more so- 
phisticated VOS models can be constructed which in- 
clude additional parameters. For example, supercon- 
ducting strings |44| carry currents described by a charge 
parameter [45], while cosmic superstrings [T2] [13] can 
be modelled as a collection of several network compo- 
nents, each with a different correlation length, velocity 
and loop production efficiency, as well as several extra 
interaction coefficients describing the formation of string 
junctions l(i IfTj. The generalisation of our methods to 
such networks will be the subject of a follow-up publi- 
cation. 

Here, we will consider one extra parameter, a, de- 
scribing the effect of short-scale wiggles on the string 
energy and tension, thus affecting the UETC. To a first 
approximation, the VOS equations do not depend di- 
rectly on a, as the maineffects of short-scale-structure 
are already captured by k in equation ^ . More sophis- 
ticated "wiggly" models exist [50], but the VOS equa- 
tions above are sufficient to reproduce the right phe- 
nomenology with respect to string correlation lengths 
and velocities. The parameter a can be readily obtained 
by Nambu-Goto simulations. Note that, physically, a 
and k must be implicitly related, but a quantitative for- 
mula is not known. 

The easiest way to describe the string wiggliness a is 
by writing the string energy-momentum tensor as 



1 For a random walk, these two length-scales effectively coincide 
and in simulations of the simplest models of strings they are 



found to be in good numerical agreement. 
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(9) 



where U is the string energy per unit length and T the 
string tension. Lorentz invariance requires U = T = /i. 
However, if short scale wiggles cannot be resolved, this 
has the effect of increasing the string energy per unit 
length and reducing the tension [5T1I52] in such a way 
that their product remains constanir 

UT = [j? . (10) 

Then, an effective "coarse-grained" energy-momentum 
tensor can be written in the form of ([9]) with "renor- 
malised" string energy density and tension, parame- 
terised through a as follows 



U = afi 



(11) 



As we will see, this parameter a enters directly into the 
UETC through this effective energy-momentum tensor 
for "coarse-grained" string segments. 



B. Energy- momentum tensor from a single string 

Let us consider a single straight segment of the string 
network. We are free to orientate the wave vector as 
k = kk 3 , where k is its magnitude, in which case the 
real component of the Fourier transform of the energy- 



momentum tensor ^ becomes (e.g 
fxa sin(kX 3 ^T /2 



) 
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e, 



v 2 X i X j 



kX 3 /2 



cos(x + kX 3 vT^ ,(12) 
-XjXi 0qo , (13) 



where X and X are randomly orientated unit vectors 

satisfying X ■ X = and i,j = 1 ... 3. The phase of 
the Fourier mode is set by the location of the string, 
Xo, where x = k ' x o- One can then identify the scalar, 
vector and tensor anisotropic stress by 

s = (20 33 -0 n - 22 )/2, 

© V = 0i3, 
T = 12 . 

For each string segment this gives 



20 s 
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0^ 

0oo 
0^ 

©00 



3X,X, - 1 - 



1 - v 2 



v 2 X ± X 3 
v 2 X l X 2 



1-v 2 



-XiX 3 



-XIX2 



(14) 
(15) 
(16) 



:5V,.V, -1) (17) 
(18) 
(19) 



r 



We choose X to be the usual position vector in spheri- 
cal coordinates, i.e. X T — (sin cos 0, sin ^ sin 0, cos 6*). 

Given the orthogonality of X and X, once you specify 

X there is an angle ip from to 2tt which gives X, 



X 



cos 6 cos 4> cos ip — sin <f) sin ip 
cos 9 sin 4> cos ip + cos 4> sin ip 
— sin 9 cos ip 



(20) 



We can then generalise to a network of strings com- 
prising many segments with different orientations, by 
averaging over angles. 



C. Unconnected segment model 

The Unconnected Segment Model (USM) frame- 
work [25] [26] models the string network as a set of uncor- 
rected straight segments, each moving with random ve- 
locity. All segments are assumed to be produced at some 
fixed initial time. Throughout cosmic history a certain 
fraction of these segments decay at each epoch in order 
to maintain scaling of the network. Since, during scal- 
ing, the number density of strings falls as n(r) oc r~ 3 , 
one needs to track an extremely large number of seg- 
ments in order to have of order one segment remaining 
today. 

To avoid tracking each segment, the USM consolidates 
all string segments that decay at the same discretised 
epoch Tj into a single string. Specifically, the number of 
segments which decay between conformal time Tj_i and 
n is 



N d (T i )=V[n( n - 1 )-n(n)) 



(21) 



where V is the simulation volume and n(r) the number 
density of strings, given by 



n(r) 



C(r) 
(£r)3 



(22) 



The coefficient C(r) is determined by requiring the to- 
tal number of strings at any time be equal to V/(£,t) 3 . 
The Nd(Ti) string segments are then consolidated into a 
single segment, which has a energy- momentum weight 
•v/ iVd(Tj). The total energy-momentum tensor is a sum 
of these consolidated segments, given by 



A" 



M „ = £ [NMt'^i^^ruLf) , 



(23) 



2 This is also the case if Lorentz invariance is broken at a fun- 
damental level, e.g. in superconducting strings |53| . 



where Q' is the energy-momentum tensor of a sin- 
gle segment i, K the number of consolidated segments 
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and T (r, r i; Lf) is a function that controls the rate of 
string decay. This has the form 



1 T<LfTi, 
r""(r.r,,I,)= \ §+l( y 3_ 3y ) L fTi <T<Ti, 

T > Tj , 

(24) 



where 



v _ MLfn/T) 



(25) 



Segments start to decay at LfTi and have disappeared 



completely at r^. The parameter Lf < 1 therefore con- 
trols how fast the segments decay, and approximates 
the Heaviside step function when Lf — > 1. The value 
of this parameter is less well understood from Nambu- 
Goto simulations. In general, higher Lf results in less 
power in the CMB power spectrum US] for a fixed 
string tension Gp,. 



D. Analytic expressions 

To compute the unequal time correlator (UETC) an- 
alytically we integrate over all strings in the network 



(e(k, Tl )Q(k,T 2 )) = 



2/(r 1 ,r 2 ,e,L / ) 
16tt 3 



2tt 



p2ir r2ix 

sinHafi I dt/j / dxe(fc,Ti)6(fc,T 2 ) . 



(26) 



where by 9 we mean ©oo, @ s , © y , © T and 
fi T iT T 2,£,,Lf) is a scaling factor associated with 
the decrease in the number density and the decay of 
strings (see below). The factor of two in the numerator 
arises from us only considering the real part of the 
UETC. Three of the integrals can be performed to give 
the remaining 9 integral in terms of Bessel functions. 
Some of the resulting terms do not have compact 
analytic expressions, so we need to make use of series 
expansions. In each case we can write the 9 integral as 
the sum 

(9(fc,r0e(fe,r 2 )) = /(T ^^y x (27) 



limit, 



A [h ,P)-h 0+ > P)] 



where the relevant integral identities Ii{x±, p) are given 
in Appendix [Aj p = k\r\ — T2\v, — k^Ti^ and 

x± = (xi ± x^)/2. For each scalar, vector and tensor 
UETC we write down the amplitude Ai of each integral 
component, in Appendix |A") 

We will also be interested in certain asymptotic limits 
of the UETC. For completeness we write down the small 



<e(fe,r 1 )e(fc,r 2 )) « k 2( 1 _ v 2) B - 



(28) 



Finally, we are also interested in the equal time correla- 
tor (ETC), so that x\ — x 2 = x + = x and £_ = p = 0. 
In this case we write the exact expression 



(9( t , T )e (t .T)) = ^;^f c. 



(29) 



The B and C coefficients are also given in Appendix [5] 
Equations ( 27 1 , ( 28 1 and ( 29 1 form the basis of our fast 



UETC code. 



E. Scaling factor 



In the above we introduced the scaling factor 
f( T i, t~2,£,, Lf). If we assume the network consists of 
only a single string whose number density per simula- 
tion v olu me is fixed and the string does not decay, then 
from (23) in the USM 6 M „ = ©L,, while analytically we 
can set j{t\, t 2 , £, Lf) = 1. Next, we turn on the scaling 
and decay of strings. This means in the USM the UETC 
is given by 



A" 



(6^(fc, Tl )Q pa {k, r 2 )> = (@l u (k, n)e^(*, r 2 )) ]T N d {n)T° s {T 1: n,Lf)T oS {T 2 ,n,Lf) . 



(30) 



I 

since (G 1 © 3 ) is equal to (0 1 © 1 ) for i = j and for i ^ j. scaling factor 
Increasing the number of segments such that K — > oo 

then the sum when Lf — > 1 can be evaluated to give the f(Ti,T 2 ,^, Lf — !• 1) 



[eM««(Ti,7s)r 



(31) 
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The < Lf < 1 scaling factor is much more lengthy 
but is possible to write down analytically. For the pur- 
poses of our work we assume the strings segments decay 
instantaneously and set Lf — > 1. 



F. Comparison with simulations 



The expressions for the UETC's given in Appendix [X] 
have been coded in a self-contained Fortran 90 module. 
In this code, we make use of asymptotic limits to im- 
prove speed and accuracy. For small x\ and X2, we use 
the small x expansion (the B coefficients) , thereby elim- 
inating the need to evaluate trigonometric functions or 
perform a series expansion of spherical Bcsscl functions. 
For x\ ~ X2, we use the form of the ETC (the C co- 
efficients), as the amplitudes Aj contain terms oc l/p 2 , 
which diverge when x\ = X2 (this divergence in the A4 
coefficients cancels when fully expanding out the inte- 
gral in this limit). 

In all other cases we perform the sum of integral com- 
ponents, tuning the number of terms in the spherical 
Bessel series expansions which are required for sufficient 
accuracy, given x\ and X2 (in general larger x\ and X2 
require more terms) . A slight caveat here is that the in- 
dividual terms of the series expansion can become much 
larger than the total sum, and for x\ or X2 greater than 
~ 30, loss of accuracy can occur due to the limitation of 
double precision arithmetic in our code. In these cases 



we perform the integrals of (Al) and (A4) numerically. 



The result is a code which can generate the scalar, vector 
and tensor UETC's, over all scales of interest relevant 
for the CMB, on a single CPU in - 20-30 seconds. We 
have parallelised this code, so a moderately threaded 
CPU can calculate the UETC's extremely quickly. 

For our comparisons with the (numerical) USM we use 
the publicly available CMBACT code [26], using a string 
decay parameter of Lf = 0.99 and K — 500 consolidated 
segments. We assume this value of Lf is close enough 



to unity so Eqn. (31 ) is valid in the analytic calculation. 
We use string network parameters v = 0.65, £ = 0.13 
and a = 1.9, turning off the time evolution of these in 
CMBACT. These parameters are chosen to match the VOS 
evolution in the radiation era, as in Refs. [26j [54] . At 
present, we have not attempted to incorporate time de- 
pendence of the string parameters in our code, as this 
complicates the eigen-decomposition of the UETC (see 
below) and obscures some of the physics of our method. 
Although the scaling solution will vary somewhat as the 
Universe changes from radiation to matter domination, 
for the purposes of our model we fix the network param- 
eters to be constant. 

In Fig. [T] we show scaled UETC's, defined as 
{t\T2) 1 ^ 2 (0(fc, Ti)0(fc, T2))//! 2 , with the numerical and 
analytical results shown respectively in the left and right 
columns. From our analytic expressions for the UETC 
one can see that the only dynamical variable is the com- 
bination fcr, which is in agreement with physical expec- 
tations. We plot the UETC derived from both 20000 




i°g 10 |ki,) 




/ 






FIG. 1. Scaled unequal time correlators (UETC's), 

(tit 2 ) 1/2 (e(fc,r 1 )e(fc,r 2 ))/ A t 2 . (Left) The UETC is derived 
from 20000 realisations of the string network using CMBACT, 
with a decay parameter Lf = 0.99, and (right) the analytic 
calculation (see text). From top-to-bottom we show the 00 
component, then scalar, vector and tensor anisotropic stress 
components, and finally the 00 — 5* cross correlator. 



realisations of CMBACT and our analytic expressions, and 
find the two in excellent agreement. For comparison, 
deriving the UETC from 20000 realisations of CMBACT 
takes several hourf0 



3 This does not include the Einstein-Boltzmann solver part of 
CMBACT, which takes much longer to run. 
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III. CMB ANISOTROPIES 

An important source of difficulty in the computation 
of CMB power spectra from cosmic strings is that they 
are incoherent sources of perturbations. As a result 
one faces complications in calculating the ensemble av- 
erage of the CMB transfer functions. In general there 
are two approaches to address the issue of incoherence. 
The first avoids dealing with the UETC directly. In- 
stead, one creates an ensemble of source histories of the 
USM, which has the same 2-point correlation statistics 
as the UETC. The source histories are then used in an 
Einstcin-Boltzmann code, and the power spectrum can 
be obtained by averaging over many realisations of the 
network. This is the approach CMBACT uses. 

The second works directly with the UETC, which 
can be estimated from string simulations, or in our 
case analytically. The UETC is first decomposed into 
eigenmodes. Each of the individual modes is coherent, 
and can be used as source functions in the Einstein- 
Boltzmann code. This approach was originally proposed 
in Ref. jSHj and is the method used by groups calculat- 
ing CMB spectra from Abelian-Higgs string simulations 
(e.g. [SI])- Since we readily have the UETC, we follow 
the latter approach. 



along the diagonal, and 00 — S as the off-diagonal com- 
ponents. After performing the eigen-decomposition the 
first half of the eigenvector corresponds to density, and 
the other half the anisotropic stress. 

This diagonalisation procedure introduces a change 
of basis. Since the eigenvectors are orthogonal, and the 
Einstein-Boltzmann solver uses linear perturbation the- 
ory, one can use each of these new basis functions as 
a source term. This means we can simply modify the 
Einstein equation sources according to 



e(fcr) 



u(kr) 



(33) 



The total power spectrum is then found by summing 
over all eigenmodes 



Ci — XiC\ , 



(34) 



where C\ are the individual power spectra of each eigen- 
mode. In practice we order modes from the highest to 
lowest eigenvalues and truncate the sum at a finite num- 
ber of modes (which gives the desired accuracy for Cg). 
Although each of the eigenmodes is coherent, good con- 
vergence can be obtained including only relatively few 
modes in the sum. 



A. Eigenmode decomposition 

Since the scaled UETC is only a function of fcr we 
can discretise it on a logarithmic grid in fcr with n x n 
grid points. It is also real and symmetric and hence 
diagonalisable. We therefore decompose it as a sum of 
eigenvectors 

(k 2 T lT2 y (T!T 2 ) 1/2 (e(fc,r 1 )8(fc,r 2 )) = (32) 

n 

^ XiUi(kri) (g) Ui(kT 2 ) , 

i=l 



where Xi is the eigenvalue and eigenvectors Ui are also 
only a function of kr. The term (k 2 TiT2) 1 is a weight- 
ing factor, designed to improve the accuracy of the re- 
constructed UETC when the eigenmode sum is trun- 
cated at some value less than n (see below). A weight 
of 7 > 0, for example, will improve the reconstruction of 
the kr > 1 region. In practice this value is tuned based 
on the scales which contribute to CMB anisotropies. We 
find a weight of 7 = 0.25 is around optimal for improv- 
ing the convergence of the power spectra. 

Since the correlation between scalars, vectors and ten- 
sors vanishes we can perform the diagonalisation on 
the vector/tensor UETC's independently. However, for 
the scalars we have two components of the energy- 
momentum tensor, the density, 0oo, and anisotropic 
stress, , which are correlated. In order to perform 
the decomposition for scalars we therefore discretise the 
UETC on a In x 2n grid. This has a 2 x 2 block diagonal 
form, with the 00 — 00 and S — S components of (32) 



B. CMB power spectrum 

The Einstein-Boltzmann solver we use is the publicly 
available CAMB code [56]. This has several advantages 
over the CMBFAST [57] code which CMBACT is based on. It 
has been tuned for efficiency and accuracy, and also uses 
the OpenMP parallel framework. Our implementation 
is based on the January 2012 version of the code. 

The string eigenvector sources were next incorpo- 
rated. For scalars we also require the longitudinal com- 
ponent of the velocity perturbation, Qrj- This can be 
derived from the known components by solving the dif- 
ferential equation [2"5] 



-me. 



3 



Op -Oo 
H 



e 



26* 



, (35) 



where overdots now represent differentiation with re- 
spect to conformal time and % = a /a. 

The CMB power spectra do not require any initial 
spectrum to be specified, since the sources are active, 
and are given by 



7T 



k 2 dkK i {I) (k,T )A 1 e {I) (k 1 To), 



(36) 



where I = S,V,T for the scalar, vector and ten- 
sor contributions, and A^(fc,To) are the tempera- 
ture/polarisation transfer functions for each mode. 
These transfer functions can be found in detail, for ex- 
ample, in Ref. [55] . 
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FIG. 2. Comparison with CMBACT. On the top panel we show scalars, in the middle vectors and on the bottom tensors. From 
left-to-right we show the TT, EE, BB and TE power spectra. Results from CMBACT using 2000 network realisations are shown 
in solid black, and from our code using the first 16 (dotted red), 32 (dashed blue), 64 (dot-dash green), 128 (dot-dot-dash 
magenta) and 256 eigenmodes (long-dashed cyan). We use string parameters of G/i = 2 x 10 -7 , v = 0.65, £ = 0.13 and 
a = 1.9. 



For our analysis we use the mean WMAP 7-year cos- 
mological parameters [S3] - that is a baryon density 
Q b h 2 = 0.02249, cold dark matter density VL b h 2 = 
0.1120, Hubble parameter H = 70.4kms _1 MpcT 1 , 
and an optical depth to reionisation r = 0.088. We 
switch off lensing of the CMB in order to directly com- 
pare results with CMBACT. The same "default" set of 
string parameters of v = 0.65, £ = 0.13 and a = 1.9 
are used, again turning off the time evolution of these 
and setting the decay parameter to Lf = 0.99 in CMBACT. 
We fix the string tension to G/i = 2 x 10 -7 , a limit which 
is allowed, for the default parameters, by current CMB 
data [2"T] . 

As in the CMBACT implementation, we found it nec- 
essary to adjust internal accuracy parameters of CAMB 
to achieve the desired accuracy in the power spectrum. 
We varied the number of terms in the Boltzmann hi- 
erarchy and the t mode spacing to compute the spec- 
trum, but the only adjustment we found necessary was 
to increase the number of k modes for the CMB source 
functions. For the aficionados sufficient accuracy was 
achieved by setting high_accuracy_def ault = T and 
accuracy Jboost = 1.5. 



The TT, EE, BB and TE power spectra are shown in 
Fig. [2j computing the UETC's on a 1024 x 1024 grid in 
the range kr — [10~ 4 , 10 4 ], with logarithmic grid spac- 
ing, and a weighting factor of 7 = 0.25. The range and 
grid size of the UETC were chosen so that the spectrum 
is insensitive to increasing the range or number of grid 
points. The weighting factor has been chosen to improve 
the convergence of the spectrum for a finite number of 
modes. 

We plot spectra truncating the eigenmode sum at 16, 
32, 64, 128 and 256 modes, and also show the results 
from CMBACT using 2000 network realisations. The con- 
vergence rate depends on the particular spectrum of in- 
terest. In general, polarisation spectra converge more 
quickly since (aside from the reionisation induced sig- 
nal) they are sourced only at last scattering. This means 
that a smaller range of scales (and hence eigenmodes) 
contribute to the source function, and the weighting fac- 
tor is chosen to preferentially reconstruct these modes. 
The vector B-mode spectrum, for example, converges in 
as little as 32 eigenmodes. The scalar temperature sig- 
nal, on the other hand, is sourced both at last scattering 
and along the line-of-sight through the Integrated Sachs 
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FIG. 3. String parameter dependence of the total TT (top) and BB (bottom) spectra. Results are derived using 256 
eigenmodes and a weighting of 7 = 0.25. In the left panels we set £ = 0.13, a = 1.9 and vary v as 0.2 (solid black), 0.4 
(dotted red), 0.6 (dashed blue), 0.8 (dot-dash green). In the middle we set £ = 0.13, v = 0.65 and vary a as 1.0 (solid black), 
1.5 (dotted red), 2.0 (dashed blue), 2.5 (dot-dash green). On the right we set v — 0.65, a = 1.9 and vary f as 0.1 (solid 
black), 0.2 (dotted red), 0.3 (dashed blue), 0.4 (dot-dash green). In each case Gfi — 2 x 10 -7 . 



Wolfe effect. This means a wider range of eigenmodes 
contribute to the spectrum. 

Plotting the difference of each spectra we find that 1% 
level accuracy in the temperature power spectrum can 
be achieved with 256 modes. On a 3.1 GHz Intel Xeon 
CPU with 8 threads, our code runs in ~ 15 minutes. 
In comparison, around 2000 network realisations are re- 
quired for CMBACT to achieve the same accuracy - on the 
same CPU, running with only the single thread that is 
possible, the computation time is ~ 30 hours. If accu- 
racy is only required in the temperature power spectra 
at the 10% level, such that 64 eigenmodes are required, 
our code takes ~ 2.5 minutes with 16 threads. Given 
that strings are constrained to < 10% of the total tem- 
perature anisotropy of the CMB, this level of accuracy is 
tolerable for parameter estimation studies. For the first 
time it may be possible to explore the string param- 
eter space using Markov-Chain-Monte-Carlo (MCMC) 
methods. 

Generally, the results from our code and CMBACT agree 
very well, considering the two different implementations. 
The only differences of note are a slight excess in power 
in the CMBACT scalar and tensor temperature fluctua- 
tions near the peak of the spectra, and a slight excess 
in the low I tensor E-mode polarisation. We have been 
unable to reconcile these small differences by changing 
details of the UETC grid, accuracy or increasing the 
number of modes. 



C. Parameter dependence 

With our code in hand, we can now begin to explore 
the string model space. Some of these parameter depen- 
dences have also been considered in Refs. [231 |2"51 150] . 
The default model is specified by v = 0.65, £ = 0.13 and 
a = 1.9. In each case we fix two of the parameters while 
varying the other. We again compute the UETC's on 
a 1024 x 1024 grid in the range kr = [10~ 4 , 10 4 ] and a 
weighting of 7 = 0.25, truncating the sum at 256 modes. 

In Fig. [3] we show the total (i.e. scalar + vector + 
tensor) TT and BB power, normalising the value of Gfi 
such that the TT power at I = 10 is the same as the 
default model. Given the default model has a tension 
G/J.Q = 2 x 10~ 7 , the ratio of \i required for the same 
C w is m/mo = [1.99,1.57,0.11,0.65], for v = 0.2, 0.4, 
0.6 and 0.8 respectively. For a the ratio is fi/fio = 
[1.50,1.20,0.95,0.78] for a = 1.0, 1.5, 2.0 and 2.5. Fi- 
nally, for £ the ratio is fi/n = [0.84,0.33, 1.79,2.25] for 
£ = 0.1, 0.2, 0.3 and 0.4. One can see that G/i has a 
strong dependence on the values of the string parame- 
ters. 

We leave a detailed exposition of the dependence of 
the power spectra on the string parameters to future 
work, but some of the broad features can be understood 
by considering the analytic form of the UETC given 
in Appendix [A] The low I amplitude is determined by 
the small x limit (the B coefficients). For scalars, this 
is complicated by the fact that density, velocity and 
anisotropic stress all act as sources. For vectors and 
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tensors, however, the low £ dependence of both temper- 
ature and polarisation (ignoring the reionisation induced 
polarisation signal) scales according to 

c i « ^) ^ 2 (l-« 2 ) • 

(37) 

As far as we are aware this explicit parameter depen- 
dence has not been written down before. We have tested 
this with the actual C/s from CAMB, and find the agree- 
ment for polarisation from £ ~ 40 to several hundred 
(where the spectrum begins to turn over) to be ex- 
tremely accurate. The relationship is less accurate for 
temperature (within a factor of two for the range of pa- 
rameters we consider), since it is also sourced by smaller 
scales along the line-of-sight. 

Other features of the spectrum can also begin to 
be understood using the analytic form of the UETC. 
Changing the velocity, for example, has a larger affect 
on the B-mode peak position than temperature. The 
scalar contribution to the latter is dominated by the 
density fluctuations near the position of the peak. Ex- 
amining the UETC coefficients, one can see that velocity 
has no effect on the amplitude of these coefficients. In 
particular, along the diagonal the velocity has no effect 
on the scale at which the UETC turns over, as can be 
seen from the form of the ETC (the C coefficient). The 
scale at which the vector anisotropic stress UETC turns 
over, however, is effected by the velocity. Expanding 
out the coefficients one can see that a higher velocity 
causes the UETC to turn over on larger scales, which is 
consistent with the shift in the B-mode peak position. 

Note that this exploration of parameter dependencies 
by varying one network parameter at a time is strictly 
not allowed, as the network parameters are related 
within a given scaling solution, as given by Eqns. ([3j 
§[6|. Here we have used this as a phenomenological 
tool for singling out the dependencies of the UETC on 
the corresponding physical variables, which can help de- 
velop an intuition of the physics involved [61] . 



is computational overhead associated with creat- 
ing an ensemble of source histories. 

2. Once the UETC's have been calculated they can 
be diagonalised, with the eigenvectors acting as 
sources in an Einstein-Boltzmann code. A much 
smaller number of eigenmodes, as opposed to av- 
eraging over a number of source histories, are re- 
quired to achieve the same accuracy in the CMB 
power spectrum. 

3. Our implementation uses the CAMB Einstein- 
Boltzmann code, which has several advantages 
over CMBFAST, which CMBACT is based on. CAMB 
has been optimised for efficiency and accuracy, is 
more modular, and supports the OpenMP frame- 
work. Secondly, CAMB supports the latest features 
such as CMB lensing and bispectrum estimation, 
and is readily compatible with the Markov-Chain- 
Monte-Carlo (MCMC) package CosmoMC [62] . 

The resulting code can compute spectra, whose temper- 
ature power spectra are accurate at the 10% level, in 
only a few minutes on a multi-threaded CPU. 

The dramatic improvement in computational time 
achieved with our approach now allows us to perform 
comprehensive scans over the string network parame- 
ter space. In particular, the compatibility of our CAMB 
implementation with CosmoMC allows us to perform cos- 
mological parameter fitting for infLation+strings scenar- 
ios, marginalising over string network parameters. This 
will be the subject of a forthcoming publication. Our 
methods can also be generalised to more complicated 
networks of cosmic superstrings, comprising string seg- 
ments with different tensions that join together in Y- 
shaped junctions [61]. This will allow us to efficiently 
probe a wide range of well-motivated inflationary mod- 
els in string theory, constraining fundamental parame- 
ters [IB] , 
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IV. CONCLUSIONS 

The main achievement of this paper has been to derive 
analytic expressions (see Appendix |A| for the UETC's 
of cosmic strings in the USM. These expressions agree 
extremely well with simulating an ensemble of source 
histories. Our formalism makes the physical connection 
with CMB anisotropies more transparent, as UETC's 
can be obtained very efficiently for different assumptions 
regarding the properties of the string network. From 
this we have written a new code to compute CMB power 
spectra. We believe that our implementation improves 
that of the existing CMBACT code in three main ways: 

1. The UETC's are calculated analytically. While 
CMBACT does not directly compute UETC's, there 
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Appendix A: Integral identities and amplitudes 

This appendix details integral identities and ampli- 
tudes used in deriving the UETC's. Note there is ac- 
tually a divergent component in both the Ii(x,p) and 
Inix, p) integrals, but these cancel as we always take the 
combination Ii(x—,p) — I%(x+,p). Therefore in the fol- 
lowing relations the divergent component is discarded. 
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where j n (x) is the spherical Bessel function. 
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We now give amplitudes of the integral components 
and the B and C coefficients for each UETC. In the 
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